UNCLASSIFIED 


r 


HeoAoditced 
luf  the 


ARMED  SERVICES  TECHNICAL  INFORMATON  AGENCY 
ARLINGTON  HALL  STATId 
ARLINGTON  12,  VIRGINIA 


UNCLASSIFIED 


NOTICE:  When  government  or  other  drawings,  speci¬ 
fications  or  other  data  are  used  for  any  purpose 
other  them  in  connection  with  a  definitely  related 
government  procurement  operation,  the  U.  S. 
Government  thereby  Incurs  no  responsibility,  nor  any 
obligation  whatsoever;  and  the  fact  that  the  Govern¬ 
ment  may  have  formulated,  furnished,  or  in  any  way 
supplied  the  said  drawings,  specifications,  or  other 
data  is  not  to  be  regarded  by  implication  or  other¬ 
wise  as  in  any  manner  licensing  the  holder  or  any 
other  person  or  corporation,  or  conveying  any  ri^ts 
or  permission  to  manufacture,  use  or  sell  any 
patented  invention  that  may  in  any  way  be  related 
thereto. 


AFCRL-62-1089 


^^e  et 

^  -I-C  YORK  UNIVERSITY 

> 

S  Physics  Dept. 


4^#^! 


2E 

i — -  b . 

^  ■.,.? 


SEARCH  REPORT  NO.  CX-59 


s 

c,'".’  ' 

C:.' 

j  WJW«N»* 

’br  «dr 


Model  for  the  Statistical  Atom 
With  Non -Vanishing  Angular  Momentum 


JUDAH  L.  SCHWARTZ  and  SIDNEY  BOROWITZ 


Contract  Nos. 


AF  19(604)4555 
DA  30-069-ORD-2581 
Nonr-285(49),  NR  012-109 


DECEMBER,  1962 


NEW  YORK  UNIVERSITY 
Physics  Department 


Research  Report  No.  CX-59 
Document  No.  APCRL-62-IO89 


A  MODEL  FOR  THE  STATISTICAL  ATOM 
WITH  NON-VANISHING  ANGULAR  MOMENTUM 


Judah  L.  Schwartz  and  Sidney  Borowltz 


The  research  reported  In  this  paper  has 
been  Jointly  sponsored  by  the  Geophysics 
Research  Directorate  of  the  Air  Force 
Cambridge  Research  Center,  Air  Research 
and  Development  Command,  under  Contract 
No.  AF  19(604)4555,  Project  No.  7635, 
Task  No.  76361,  and  the  Office  of 
Ordnance  Research  under  Contract  No.  DA- 
3O-O69-ORD-258I,  Project  No.  2360,  and 
the  Office  of  Naval  Research  under  Con¬ 
tract  Nonr-285(49),  NR  012-109. 


DECEMBER  I962 


Requests  for  additional  copies  by  Agencies  of  the  Department  of 
Defense,  their  contractors,  and  other  Government  agencies  should 
be  directed  to  the: 

ARMED  SERVICES  TECHNICAL  INFORMATION  AGENCY 
ARLINGTON  HALL  STATION 
ARLINGTON  12,  VIRGINIA 

Department  of  Defense  contractors  must  be  established  for  ASTIA 
services  or  have  their  ' need- to -know •  certified  by  the  cognizant 
military  agency  of  their  project  or  contract.  All  other  persons 
and  organizations  should  apply  to  the: 

U.S.  DEPARTMENT  OP  COMMERCE 
OFFICE  OP  TECHNICAL  SERVICES 
WASHINGTON  25,  D.C. 


ABSTRACT 


In  the  first  of  the  coupled  integro -differential  equations  satisfied  by 
the  Green's  functions  of  the  many  body  system  the  Hartree-Fock  approxima¬ 
tion  is  made  for  the  two  particle  Green's  function.  The  resulting  equation  is 
written  in  a  mixed  position -angular  momentum  representation  in  such  a  way 
that  use  may  be  made  of  whatever  empirical  information  about  the  angular 
momentum  of  the  system  is  available. 

In  the  mixed  position-angular  momentum  representation  the  one 
particle  Green's  function  appears  as  a  sum  of  terms  each  of  which  corres¬ 
ponds  to  a  different  value  of  the  angular  momentum.  For  the  principal  groups 
of  the  periodic  table  in  which  the  angular  momentum  of  the  atom  is  carried 
by  either  one,,  two  or  three  electrons  added  to  or  missing  from  a  closed  shell 
the  appropriate  terms  in  the  expression  for  the  one  particle  Green's  function 
are  identified.  The  first  of  the  coupled  integro -differential  equations  for  the 
Green's  functions  is  then  solved  to  lowest  order  in  h  in  the  manner  indicated 
by  Baraff  and  Borowits.  The  result  is  a  modified  Thomas -Fermi  model  that 
differs  from  the  previous  results  of  Sessler  and  Foley  for  the  same  problem. 

The  Sessler-Foley  modification  of  the  Thomas -Fermi  atom  is  re- 
derived  in  the  spirit  of  the  cranking  model  for  the  nucleus  by  spinning  the 
potential.  The  present  formalism  is  seen  to  lead  to  the  Sessler-Foley  re¬ 
sult  in  the  limit  of  many  particles  carrying  the  angular  momentum.  How¬ 
ever,  the  region  of  validity  of  the  cranking  model  in  the  nucleus  lies  far  re¬ 
moved  from  the  magic  numbers  where  there  are  many  nucleons  outside  a 
spherical  core.  Thus  the  Sessler-Foley  result  may  be  interpreted  as  a  kind 
of  cranking  model  of  the  atom. 

In  contrast  to  the  previous  work  the  present  formalism  yields  a  sign 
change  for  the  quadrupole  coupling  constant  on  passing  through  closed  shells 
and  calculated  values  for  the  quadrupole  coupling  constant  of  the  halogens  are 
seen  to  be  within  an  order  of  magnitude  of  the  accepted  values. 
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ERRATA 

Eq.  II. 8  Jh  =  /  (rxD(r))n(r)<?r 
Line  8  Read 

Eq.  III. 7  Read  'i'pjin,(r)  =  ( 2v ) 3/2£i/p/hr  ® ^ 

Eq.  III.14  Read  for  right  hand  side  (2Trn)‘*G^(r^t^r^t2) 

The  second  line  of  Eq.  III. 19  should  read 

IpI .‘.expi/!i(p.rj_  +  a)(t^-tj^))  dr_^d(t^-t|) 

Eq.  111,20  Read  KG  instead  of  KG 

F  P 

K.B,  The  quantity  p  refers  to  the  magnitude  of  the  momentum 
unless  otherwise  indicated.  The  argument  of  the  spherical 

harmonic  is  (fi  ). 

V 

Line  11  Read  ....This,  however,  is  true  for  all  i  and  ... 

Eq.  III. 22  f’or  read  17. 

Eq.  III.2U  L  =  N(17)^'^^ 

^  max 

Eq.  III. 26  Read  6  instead  of  G 

000  ooo 

Eq.  III. 27  Read  instead  of  . 

Line  1^4  Read  ....Thus  has  the  significanre  .... 

Eq,  III. 33  The  operator  in  the  second  integral  should  be  2L(1)*L(2) 

Line  l6  Read  ...  table  have  three  p  electrons... 

Eq.  III. 45a  The  denominator  should  read  )  .  (R)  . 

xi  Lrr  pXy  - 

Eq.  111,49a  The  right  hand  side  should  read 

•  ,  yrT.r.(Ji) 
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I.  INTRODUCTION 

The  method  of  Thomas^  and  Fermi ^  allows  one  to  calculate  a  charge 
distribution  for  an  atom  having  no  net  angular  momentum.  This  is  usually 
done  by  assuming  that  the  electrons  in  the  atom  constitute  a  degenerate  Fermi 
gas  at  zero  temperature,  thereby  allowing  one  to  obtain  a  relationship  between 
the  Fermi  momentum  and  the  spatial  density.  The  potential  energy  is  expressed 
in  terms  of  the  electron  density  and  finally  the  total  energy  is  varied 
subject  to  the  condition  that  the  number  of  electrons  be  fixed.  The  result 
of  this  variation  in  conjunction  with  Poisson’s  equation  yields  a  second 
order  non-linear  differential  equation  for  the  self-consistent  potential  in 
the  atom.  The  electron  density  may  then  be  recovered  from  the  solution  for 
the  potential. 

The  Thomas-Fermi  equarion  as  ordinarily  derived  describes  an  atom 
with  vanishing  angular  momentum  since  the  assumed  a  priori  density  in  phase 
space  is  spherically  symmetric.  Sessler  and  Foley^have  extended  the 
Thomas-Fermi  model  to  include  non-vanishing  angular  momentum  in  the  hope  of 
obtaining  a  model  that  would  include  angular  ’’correlations**.  This  was  done 
in  order  to  deal  with  certain  atomic  properties  that  are  particularly  sensitive 
to  the  angular  distribution  of  electrons  in  the  atom,  in  particular  the 
quadrupole  coupling  constant. 

In  order  to  insert  an  angular  momentum  into  the  statistical  atom 
Sessler  and  Foley  introduce  a  departure  from  spherical  symmetry  in  phase  space 
by  displacing  the  Fermi  sphere  from  the  origin  by  a  position  dependent  amount. 
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There  therefore,  a  net  linear  momentum  at  each  point  in  the  atom  and  an 
angular  momentum  may  now  be  introduced  in  what  is  essentially  a  classical 
fashion.  Once  again  the  energy  of  the  atom  is  expressed  in  terms  of  the 
electron  density,  and  the  total  energy  is  varied  subject  to  the  constraints 
that  both  the  number  of  electrons  and  the  angular  momentum  be  constant. 

In  Section  II  we  shall  show  that  the  Sessle^f^^Toi-ey  modification 
of  the  Thomas-Fermi  model  may  be  derived  in  an  alternative  fashion,  i.e. 
by  rotating  the  nuclear  potential  while  leaving  the  Fermi  sphere  centered 
at  the  origin  in  momentum  space.  However,  this  is  precisely  the  way  in 
which  the  "cranking  model"  of  the  nucleus  first  introduced  by  Inglis^*^ 
to  explain  low  lying  nuclear  rotational  levels  is  constructed.  Indeed, 
many  of  the  features  of  the  cranking  model  have  direct  analogues  in  the 
Sessler-Foley  modification  of  the  Thomas-Fermi  model.  The  correct  value 
of  the  angular  momentum  is  obtained  through  a  suitable  adjustment  of  the 
angular  velocity  of  the  potential.  This  procedure,  however,  does  not  allow 
for  a  proper  description  of  the  symmetries  of  the  spatial  charge  distribution 
that  are  associated  with  a  quantum  mechanical  angular  momentum.  The  model, 
therefore,  is  not  able  to  yield  any  of  the  properties  of  the  atom  that 
vary  in  the  characteristic  way  of  the  periodic  table.  The  quadrupole 
coupling  constant,  for  example,  as  calculated  by  Sessler  and  Foley  is  always 
of  the  same  sign  and  has  the  same  z  dependence  independent  of  the  value  of 
the  angular  momentum. 


In  treating  a  quantum  mechanical  system  statistically,  one  sets  up 
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an  ensemble  of  systems  the  average  properties  of  which  are  taken  to  be 
reasonable  estimates  of  the  behavior  of  the  actual  system  of  interest. 

The  representative  ensemble  is  chosen  so  as  to  agree  with  our  partial 
knowledge  of  the  system  of  interest  and  is  otherwise  constructed  in 
accordance  with  some  postulate  as  to  the  a  priori  likelihood  of  different 
possibilities* 

In  the  case  of  the  statistical  atom  the  partial  information  used 
in  the  Sessler  and  Foley  work  consists  of  a  knowledge  of  the  energy  and 
the  total  angular  momentum.  All  possible  states  consistent  with  the  energy 
and  angular  momentum  of  the  system  are  assumed  equally  probable.  However ^ 
one  has  available  a  fair  amount  of  qualitative  information  as  to  the 
internal  distribution  of  angular  momentum  in  the  atom  which  could  be  used 
to  restrict  further  the  volume  in  phase  space  available  to  the  system. 

In  fact  for  the  elements  in  the  principal  groups  the  angular  momentum  is 
known  to  be  carried  by  either  onei  two  or  three  electrons  either  added  to 
or  missing  from  a  spherical  core* 

We  have  I  therefore ^  addressed  ourselves  once  again  to  the  problem 
of  attaining  an  approximate  description  of  the  atom  that  preserves  in 
addition  to  the  energy  some  other  quantum  number,  in  particular  the  value 
of  the  total  angular  momentum.  In  order  to  go  beyond  the  work  of  Sessler 
and  Foley  we  shall  make  use  of  the  fact  that  we  are  not  in  total  ignorance 
as  to  how  the  angular  momentum  is  distributed  within  the  atom.  Including 
this  information  enables  us  to  obtain  atomic  properties  that  vary  in  the 
characteristic  fashion  of  the  periodic  table. 
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We  shall  also  show,  within  the  context  of  the  present  formalism, 
that  in  the  limit  of  no  knowledge  of  the  internal  distribution  of  angular 
momentum  the  Sessler-Foley  result  remains  as  the  only  thing  one  can  say  about 
the  statistical  atom  with  angular  momentum. 

The  present  work  utilizes  the  approach  of  Baraff  and  Borowitz^  in 
applying  the  Green’s  function  techniques  of  the  many  body  problem  to  the  atom. 
In  their  work  a  mixed  position  momentum  representation  of  the  one  particle 
Green’s  function  was  expanded  in  powers  of  h  •  The  lowest  order  term  is 
seen  to  lead  to  the  Thomas-Fermi  density  and  higher  order  terms  yield  both 
the  exchange  and  quantum  corrections  to  the  density.  One  of  the  attractive 
features  of  the  method  is  that  no  a  priori  assumptions  are  made  about  the 
nature  of  the  density  in  phase  space. 

In  order  to  construct  a  formalism  within  the  framework  of  which  we 
may  derive  a  modified  Thomas-Fermi  model  with  non-vanishing  angular  momentum 
we  shall  introduce  in  Section  III  a  mixed  position-angular  momentum  repre¬ 
sentation  of  the  one  particle  Green’s  function.  This  is  accomplished 
through  the  use  of  a  mixed  Fourier-Hankel  transform.  The  transform  and  its 
inverse  are  defined  and  the  physical  significance  of  the  transform  is 
discussed. 


In  the  mixed  position-angular  momentum  representation  the  one 
particle  Green’s  function  appears  as  a  sum  of  terms  each  of  which  corresponds 
to  a  different  value  of  the  angular  momentum.  One  then  attempts  to  identify 
those  terms  that  are  significant  in  any  given  physical  situation.  In  the 
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event  that  the  angular  momentum  of  the  atom  is  due  to  one^  two  or  three 
electrons  outside  a  spherically  symmetric  shell  a  form  for  the  one  particle 
Green’s  function  is  obtained  that  yields  the  proper  quantum  mechanical 
expectation  value  of  the  angular  momentum.  A  modification  of  this  form 
allows  one  to  discuss  the  case  of  one,  two  or  three  electrons  missing  from 
a  closed  shell.  In  the  Hartree-Fock  approximation  the  first  of  the  coupled 
integro-differential  equations  satisfied  by  the  Green’s  functions  is 
solved  to  lowest  order  in  Pf  and  a  modified  Thomas-Fermi  model  is  obtained. 

On  the  other  hand,  if  many  electrons  contribute  to  the  angular 
momentum  it  is  shown  that  a  different  assumption  about  the  form  of  the  one 
particle  Green’s  function  in  the  mixed  position-angular  momentum 
representation  must  be  made.  The  resulting  equation  for  the  electron 
density  in  the  atom  in  this  case  turns  out  to  be  the  Sessler-Foley 
modification  of  the  Thomas-Fermi  model.  The  quantum  mechanical  expectation 
value  of  the  angular  momentum  vanishes  in  this  case.  This  is  to  be  expected 
since  the  Sessler-Foley  model  corresponds  to  the  uniform  rotation  of  a 
rigid  sphere.^ 

In  Section  IV  the  present  model  is  used  to  investigate  the  behavior 
of  the  quadrupole  coupling  constant.  The  sign  of  the  constant  is  shown  to 
be  given  correctly  and  the  value  of  the  constant  is  shown  to  increase  with 
z  in  the  case  of  the  halogens.  The  calculated  values  are  within  an  order  of 
magnitude  of  the  accepted  values. 


II  THE  SESSLER-FOLEY  AND  CRANKING  MODELS 


In  their  work  Sessler  and  Foley  arrive  at  the  conclusion  that  a 
statistical  atom  with  angular  momentum  rotates  as  a  whole.  This  result  i 
consistent  with  a  theorem  proven  by  Landau  and  Lifshitz®  that  states  that 
the  only  macroscopic  motions  compatible  with  an  equilibrium  state  of  a 
statistical  system  are  translation  with  a  uniform  velocity  and  rotation 
with  uniform  angular  velocity  about  a  fixed  axis.  This  can  be  shown  in 
the  following  fashion: 

Let  the  system  be  subdivided  into  small  parts.  Each  part  is  to 
be  sufficiently  large,  however,  that  an  entropy  function  may  be  defined 
for  it.  Let  mj_,  and  E^  denote  the  mass,  momentum  and  energy  of  the 
if^  part  of  the  system,  respectively.  The  entropy  of  the  system  *  then 

II.  1  s  =  Si(Ei  -  Pj^) 

i  2m£ 

In  addition,  for  a  closed  system  we  have, 


II. 

2a 

1  -  CONSTANT 

i 

II. 

2b 

1  Ti  a  CONSTANT 

i 


We  now  introduce  the  Lagrangian  multipliers  a  and  b  and  maximize 
the  entropy  subject  to  the  constraints  that  the  momentum  and  angular 
momentum  of  the  system  be  conserved. 
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The  variation  yields  the  result  that  for  maximum  entropy 
^  li  =  aT  -  (Fi  X  FT) 

mi 

Khere  T  is  the  temperature*  We  see,  therefore,  that  the  only  macroscopic 
motions  consistent  with  a  statistical  system  in  equilibrium  are  linear 
translation  with  velocity  aT  and  uniform  rotation  with  angular  velocity  bT* 

This  result  depends  on  the  fact  that  we  have  no  knowledge  of  the 
internal  distribution  of  angular  momenta;  We  know  only  its  total  value. 

One  might,  therefore,  anticipate  a  different  result  if  one  knew  on  empirical 
grounds  that  the  entire  angular  momentum  of  the  system  was  carried  by  a 
single  particle.  This,  then,  would  correspond  to  a  further  restriction 
of  the  volume  in  phase  space  available  to  the  microcanonical  ensemble.  We 
shall  indeed  find  that  one  is  led  to  different  equations  for  the  density 
in  the  system  when  one  inserts  empirical  conditions  known  to  hold.  Indeed, 
even  without  the  inserting  of  empirical  information,  a  detailed  treatment 
allowing  for  internal  degrees  of  freedom  can  be  made  leading  for  different 
results  than  those  illustrated  above, 

It  would  therefore  seem  reasonable  to  find  that  the  statistical 
atom  with  only  the  value  of  the  total  angular  momentum  specified  should 
rotate  rigidly  as  indeed  was  found  by  Sessler  and  Foley,  These  workers 
introduced  an  angular  momentum  into  the  Thomas  Fermi  atom  where  it  is 
ordinarily  absent  by  displacing  the  Fermi  sphere  in  momentum  space  from 
the  origin  by  a  position  dependent  amount  D(r),  This  imparts  a  net  linear 
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momentum  to  each  point  in  the  atom  and  a  non-zero  angular  momentum  is  now 
possible.  The  calculation  is  outlined  briefly. 


The  kinetic  energy  is  given  by 

'  (  f 

II.  4  ^1  =  2  dr  dp  p^  =  2  dr  dp  (pfD)^ 

T2WFj~J-2sr 

disp.  sphere 

=  2  dr  dp  (p2+D2+2^*D)  =1  4r  f  Pf  f 

TTirKP'  2m  J  '‘1^ 


direction I 
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One  now  introduces  the  Lagrangian  multipliers  X  and  u  on  the 
conditions  of  angular  momentum  and  particle  number.  Variation  with  respect 
to  D(r)  yields 

II 4  9  =  -Xmr  sin  0 

Variation  with  respect  to  n(r)  in  conjunction  with  Poisson’s  equation  yields 
the  following  equation  for  the  potential  in  the  atom; 


II.  10 

where 
II.  10a 


V^<J>  =  4e^  [2m(<J>~p)  +  (Xmr  sin  6 


n(r’)dr*  =  potential  at  r 

ff“£T 


The  quantity  y  corresponds  to  the  chemical  potential  and  vanishes 
in  the  case  of  the  neutral  atom. 

The  form  of  the  modified  Thomas-Fermi  equation  II.  10  indicates 
that  the  model  describes  a  rigid  sphere  rotating  with  constant  angular 
velocity  X ,  and  the  correct  value  of  the  angular  momentum  is  obtained  by 
a  suitable  adjustment  of  the  value  of  X.  ThiS|  however,  is  a  classical 
treatment  of  the  angular  momentum.  To  begin  with,  the  total  angular 
momentum  is  taken  equal  to  the  z  component  in  the  above  derivation  and 
further,  quantum  mechanically  it  is  meaningless  to  speak  about  an  angular 
momentum  for  a  spherically  symmetric  system  since  the  wave  function  of 
such  a  system  is  invariant  under  rotation. 

This  modification  of  the  Thomas-Fermi  model  may  be  derived  in  an 
alternative  fashion  that  is  very  much  in  the  spirit  of  the  cranking  model 
of  the  nucleus.  This  model  was  first  introduced  by  Inglis  in  an  attempt 
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to  explain  the  rotational  levels  that  arise  from  the  more  or  less  rigid 
rotation  of  nuclei*  The  derivation  of  the  cranking  model  hinges  on  the 
notion  that  the  nucleons  are  in  a  potential  well  that  is  externally  rotated* 
We  shall,  therefore,  attempt  to  rederive  the  Sessler-Foley  modification 
of  Thomas-Fermi  model  by  leaving  the  Fermi  sphere  centered  at  the  origin 
in  momentum  space,  and  rotating  the  nuclear  potential  well  instead* 


The  kinetic  energy  of  the  atom  is  now  given  by 


II.  11 

Ek  =  2 

dr 

dp  p2  =  1 

(2nh)3  , 

o  , '^'m  mh^  , 

Sphere  ' 

dr  Pf5 


Rotating  the  nuclear  potential  well  will  have  the  effect  of  adding 
a  centrifugal  term  to  the  potential  energy  which  becomes 


II.  12  E  = 

P  2 


dr' 

dr  n(r)n(r ' ) 

J  j 

ir'-ri 

dr  n(r)  ^eV  (r)  -  mr^w^j 
^  2 


The  energy  density  E'  is  related  to  the  energy  density  E  of  the 
non-rotating  system  by 

II.  13  E'  =  E  t  n(r)mr^ai2 

2 

We  may  therefore  consider  the  quantity 
II.  14  Cmr^ai]n(£) 

as  an  angular  momentum  density.  The  integral  of  this  density  over  the  volume 
of  the  atom  is  the  total  angular  momentum  which  is  to  be  conserved. 
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Introducing  the  Lagrangian  multipliers  X  and  p  for  the  conditions 
of  conservation  of  angular  momentum  and  particle  number  respectively  we 
obtain  as  the  expression  to  be  minimized 


II.  15 


V3 

^  (  3tt^)  [n(r)]  + 

m  lOiT^  2 


dr*  n(r)n(r*)  -  eVu(r)n(r)  +  mr^w 


r-r* 


t  X 


dr  n(r)mr^w  +  P 


dr  n(r)  . 


Variation  with  respect  to  o)  taking  into  account  the  fact  that  only 
components  of  cu  normal  to  the  polar  axis  contribute  to  the  angular  momentum 
yields 


II.  16  0)  =  -XSin  0 

Variation  with  respect  to  n(r)  in  conjunction  with  Poisson’s 
equation  yields  finally 


2 

“  ^^^2  ^2m((i)-p)  +  [xmrSin 

3it 

which  is  seen  to  be  identical  with  II.  10. 


We  see,  therefore,  that  the  external  rotation  of  the  potential 
well  leads  to  the  same  equation  for  the  density  as  that  obtained  from  the 
displacement  of  the  Fermi  sphere  in  momentum  space. 

The  rotational  levels  of  nuclei  so  well  explained  by  the  cranking 
model  do  not  appear  in  nuclei  in  the  vicinity  of  the  magic  numbers.  There 
are  for  example  no  rotational  levels  in  and  it  is  only  as  one  goes 


^  n(r) 


12 


away  from  lead  toward  the  heavier  nuclei  that  rotational  levels  arise  and 
the  theoretical  level  spacing  is  approached.  This  is  also  true  in  the  rare 
earth  region  far  removed  from  the  magic  numbers.  V/e  see  therefore  that 
a  model  based  on  the  external  rotation  of  a  potential  well  has  its 
greatest  utility  when  applied  to  situations  in  which  large  numbers  of 
particles  lie  outside  a  spherical  core.  We  shall  in  the  following  sections 
attempt  to  discover  whether  this  is  true  in  the  case  of  the  atom  as  well. 
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III  -  THE  MIXED  POSITION  -  ANGULAR  MOMENTUM  REPRESENTATION 

V/e  use  the  Green’s  function  formalism  introduced  by  Schwinger^  for 
the  quantum  mechanical  many  particle  system  in  the  ground  state.  This  des¬ 
cription  of  the  many  body  system  is  a  time  dependent  one  resembling  the 
time-independent  description  of  the  many  body  system  employing  the  set  of 
1,  2,,,,N  particle  density  matrices.  In  this  formalism  the  N  particle 
Green’s  function 


is  analogous  to  the  N  particle  density  matrix  and  contains  the  same  informa¬ 
tion  about  the  system, The  Green’s  functions  also  satisfy  a  hierarchy  of 
equations  in  which  the  N  particle  Green’s  function  is  coupled  to  the  (N+1) 
particle  and  the  (N-1)  particle  Green’s  functions.  The  procedure  to  be 
adopted  in  v;hat  follows  is  that  of  Baraff  and  Borowitz  in  which  a  suitable 
representation  of  the  one  particle  Green’s  function  is  expanded  in  powers 
of  h. 


The  hamiltonian  for  an  atom  having  a  nucleus  of  infinite  mass  and 
charge  Ze  at  the  origin  surrounded  by  N  electrons  is  given  by 
III.  la  H  =l  Ho(r^)  +  1/2  Iv(r^-r.) 

i  i/j 

III.  lb  H^(r.)  =  v?-  Ze2 

III.  Ic  v(r,--r^)  =  e2 

— J  - - ^ 


The  first  of  the  coupled  hierarchy  of  integro-differential  equations 


involving  the  one  and  two  particle  Green’s  functions  is 
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HI.  2  (iha_-  Ho(rj_))G^(r^t^:r’t|)ti/dr^v(r^-r2)G2(r^t^r2t^:r|t{r2t|} 
811  ’ 


=  h6(r^-r|)5(tj^-t|) 


where  the  coordinates  in  the  arguments  of  the  Green’s  functions  are  assumed 
to  include  spin  direction  as  well  as  spatial  position,  and  summation  over 
the  corresponding  spin  index  is  assumed  to  accompany  integration  over  a 
coordinate. 


We  introduce  the  Hartree-Fock  approximation  by  replacing  the  two 
particle  Green’s  function  by  an  antisymmetrized  product  of  one  particle 
Green’s  functions,  i,e. 


Ill,  3  ^2 :rj_tir^t^)  ,  :£iti )ei(r2t2 : ) 

-Gi(i^ti:r^tpGi(r2t2:riti) 


If  we  take  note  of  the  fact  that  in  the  absence  of  spin  dependent 
forces  a  particle  propagates  without  changing  its  spin  we  may  insert  III. 3 
in  III, 2  and  perform  the  required  summation  over  spin  indices  obtaining 

III.  4  (ihj.  -  Ho(rjL)]Gj_(rj^ti:r|tp 

ati 

+i/d£2''(£l-£2)^2Gi£r2ti:r2t+)Gi(riti:r^tp-(^(r^t^:r^tpG^(r2t^;r^tp) 


=  h6(r  -r' )6(t  -t ' ) 
-1-1  11 
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Following  Baraff  and  Borowitz  this  may  be  written  as  a  pair  of  in¬ 
tegral  equations: 

III.  5a  /^(£l'ti!r2t2)dr2dt2Gj^(r2t2r|tp  =  k6(r-£|)6(t^-t^) 

where  the  kernel  ^  defined  by 

III.  5b  ^(£I^1’£2^2^  “ 

(ih  3  -  ^o(^)  +  2i  Jd£'v(r-r’ )G2^(£*t:£*t'*’)  ]^^£]_“'£2^^^^1~^2^ 
a^i  ^ 

-iv(rjL‘*£2  )^i(£l‘^1*£2'^"*' ^ 

It  is  at  this  point  that  Baraff  and  Borowitz  effect  a  transformation 
to  a  mixed  position-momentum  representation.  For  the  description  of  an  atom 
with  non-vanishing  angular  momentum  it  will  prove  useful  to  choose  a  more 
appropriate  representation.  The  importance  of  the  choice  of  representation 
rests  on  the  fact  that  the  many  body  problem  we  have  formulated  is  not 
solved  exactly.  If  it  were,  then  all  the  properties  of  the  system,  including 
of  course  the  angular  momentum  could  be  deduced  from  the  solution  independent 
of  the  representation  of  the  Green*s  function.  However,  one  is  rarely  so 
fortunate  as  to  be  dealing  with  a  problem  that  is  soluble  exactly.  If,  then, 
as  we  have  indicated,  we  wish  to  take  the  point  of  view  that  we  are  going 
to  specify  in  our  approximate  solution  something  of  the  internal  behavior 
and/or  structure  of  the  system  that  we  may  know  from  other  sources,  we  must 
use  a  representation  that  is  appropriate  for  the  inclusion  of  such  additional 
information.  V/e  introduce  therefore  the  mixed  position  angular  momentum 
representation  which  we  have  found  to  be  most  useful  for  these  purposes. 

By  analogy  to  the  Fourier  transform  obtained  by  multiplying  the 


16 


function  to  be  transformed  by  the  solution  to  the  free  particle  wave  equation 
in  Cartesian  coordinates,  followed  by  an  appropriate  integration,  we  shall 
introduce  a  transform  that  makes  use  of  the  free  particle  solutions  in 
spherical  coordinates,  The  usefullness  of  this  procedure  is  apparent  when  one 
remembers  that  a  Fourier  transform  decomposes  a  spatial  description  of  a 
particle  into  plane  waves  each  corresponding  to  a  particular  linear  momentum# 
The  procedure  adopted  here  will  produce  a  decomposition  into  spherical  waves 
each  one  of  which  carries  a  particular  angular  momentum.  Since  the  transform 
is  somewhat • unusual  we  shall  describe  it  in  detail. 

The  single  particle  Green’s  function  is  to  be  transformed  as  follows: 


where 


III.  7  H-  (r)  =  (2Tr)3/2n  j  (pr)  ^  (g  ,) 

-  (4[r  Z+1/2  n-" 

where  cylindrical  Bessel  functions,  the  ,4)  )are  normalized 

spherical  harmonics  and 


III#  8 


R  =  "  -2 


is  to  be  held  fixed  in  evaluating  the  integral  defined  by  the  transform# 

Since  the  functions  and  have  their  time  and  energy  dependence 
related  through  a  Fourier  transform,  that  part  of  the  calculation  of  the  inverse 

A 

offers  no  difficulty.  The  inverse  of  the  angular  momentum  dependence  of  G 

X/ 

is  calculated  by  multiplying  by  the  complex  conjugate  of  y  and  integrating 

p  X  m 

over  the  magnitude  of  the  momentum  and  summing  over  all  values  of  Z  and  m# 

Thus  we  have  from  III#  6 


la) 


^,m 
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The  time  and  energy  integrations  may  be  done  immediately  yielding 
for  the  right  hand  side  of  III,  9 

III.  10  2TrHZ  /g  (r  t'r  t')H'  (r')dr  dp 

Jim  1  1  12  2  pitm  —1  pZm  —1  —1 


We  may  now  perform  the  p  integration  by  noting  that 


11 


f  J  (kx)J  (kx  )kdk  =  <5(x-x  ) 

J  m  mo  o 

independent  of  the  value  of  m. 

The  result  of  the  p  integration  in  III.  10  is 

in.  11  (2.5)“  ''l»‘®ll>'L<''l*'l>“"VV*l 

Since  the  spatial  integration  is  to  be  performed  keeping  R  constant 
there  is  only  one  independent  angle  variable  and  we  may  expand  the  spatially 
dependent  one  particle  Green's  function  as  follov/s: 

III.  12  G,  (r  t'r  t' )  =  5:  G^''(r  t'r  t'  )Y*  (  0  .({.  ) 

^-11-2  2  uv  1  112  2  uv  r  1 

Expression  III,  11  then  becomes 


III 


r  r* 
1  1 


The  spatial  integrations  may  now  be  performed  yielding 

^0  A  ^ 

r  G  ^"'(r't'r't')  Y"  (0'(^')=(2  ifR)  G  (r't'r'f) 
)[tn  1  112  2  1^1  1112  2 


I  ) 

1  sin0d6d4» 
1  1 


III,  14  (2ttR) 


4 
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where  use  has  been  made  of  III.  12  in  the  last  step. 
The  inverse  transform  is  thus  given  by 


*  “’iil  /  N 

in.  15  ^  S_^JRp„)  B  <*r"2>dp<,„ 

The  spatial  density  of  the  electrons  may  now  be  expressed  in  a 
form  that  displays  the  usefulness  of  the  mixed  position-angular  momentum 
representation.  For  the  atom  with  non-zero  angular  momentum. 


III.  I6a  „<y  ^ 

III.  16b  n  (Rp)  =  -if  G  (Row)  ^ 
m  — ^  ^  £m  —  2TTh 


Equations  III.  16  suggest  the  following  physical  interpretation;  the 

quantity  n  (R,  p)isa  density  in  a  phase  space  that  is  chara^^^  i.‘ized  by  a 
j^m  ~ 

momentum  subspace  with  one  continuous  coordinate  and  two  discrete  indices 

indicating  the  symmetry  properties  of  the  angular  momentum.  Physically  we 

shall  interpret  giving  both  the  spatial  distribution  as  well  as 

the  angular  momentum  distribution  of  electrons  in  the  atom. 

The  mixed  Fourier-Hankel  transform  of  equation  III.  5a  is  taken  with 

r  tr  ^ 

respect  to  r_  and  t^-t*  keeping  R  =  — 1  ~1  constant.  We  obtain 
1  11  2 

i^(t  -t • )  t 

III.  17  /K(r  t  r  t  )dr  dt  G  (r  t  r’t*)  ^  (r  )e  ^  ^  ^  dr  d(t  (R) 

^  -1  1“2  2  -2  2  1  -2  2-’l  1  pjlm  1  “1  1  1  8  pjim  - 

1 2 

The  quantity  ^  (r  )  may  be  expanded  in  plane  waves  as  follows 
pirn  ^1 
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ip*r 

=  M  /  da  (a  ) 

n  p  £m  p 

Inserting  III.  17  in  III.  18  and  interchanging  the  order  of 
integration  we  obtain 


(R)  =  /  dr, -t,' ) 


III.  19  H 

iP£.+i±i(tT-t’ ) 

=  -^  /  dfipY^m(fip^/‘^^-1^2£2^2)^Il2^^2^l(£2^2li^P®  ^  ‘  ^rd(^-^)=y  dflY  ajp)l^G 


v/here  KG  is  a  Fourier  transform  defined  by  the  above  integral.  This 

.  .  . 

relationship  between  the  Fourier  transform  KG^  and  the  Hankel  transform 

KG  is  quite  general,  i.e. 

Pilm 


III.  20a  kg 


p^m  ^  P  Jim 


/  dn  Y  (fir,)  KG 

'  Am  P  r 


III.  20b 


/V..  i£ 


KG  =  r  d_  Y„  («„)  KG 
P  am  °  ^  pam 


Baraff  and  Borowitz  have  shown  that  the  Fourier  transform  KG  may 
be  expanded  in  terms  of  the  individual  transforms  K  and  G  by  means  of  a  dif¬ 
ferential  operator  of  infinite  order  which  has  the  structure  of  an  infinite 
series  of  terms  each  of  which  contains  an  explicit  power  of  R.  The  zeroth 
order  term  is  just  the  product  of  the  individual  transforms.  In  this  ap¬ 
proximation  therefore  III.  19  becomes 


^-3 

III.  21  ^  't' 

O 


pJlm 


Au  p 


1 
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where  explicit  use  has  been  made  of  the  expansion  in  powers  of  h  given  by 
Baraff  and  Borowitz  for  both  the  kernel  as  well  as  the  one  particle  Green's 
function. 

At  this  point  we  would  like  to  be  able  to  make  some  statement  as 
to  the  relative  significance  of  the  terms  in  the  one  particle  Green's 
function,  in  light  of  whatever  partial  knowledge  we  may  have  about  the 
system.  For  this  purpose  we  shall  invoke  certain  statistical  arguments. 

We  regard  the  number  of  electrons  as  large  and  construct  a  microcancnical 
ensemble  to  represent  the  system.  All  possible  internal  distributions  of 
angular  momentum  consistent  with  the  total  energy  of  the  atom  are  to  be 
regarded  as  equally  likely,  A  justification  of  the  use  of  the  microcanonical 
ensemble  for  the  isolated  quantum  mechanical  system  is  discussed  in  detail 
in  Tolman. 

If  our  partial  knowledge  of  the  system  is  such  that  we  know  that 
the  angular  momentum  of  the  system  is  carried  by  a  single  particle,  then  it 
is  clear  that  this  can  occur  in  only  one  way  and  as  we  shall  see,  there  is 
an  unambiguous  form  of  the  one  particle  Green's  function  that  describes  this 
situation.  On  the  other  hand,  if  we  know  the  angular  momentum  of  the  system 
is  carried  by  many  particles  the  situation  is  not  as  clear.  This  is  because 
there  are  many  possible  internal  distributions  of  angular  momentum  that 
will  give  rise  to  the  same  total  angular  momentum.  We  are,  however,  dealing 
vfith  an  equilibrium  state  of  the  atom  since  we  are  describing  the  ground 
state  and  therefore  we  shall  adopt  the  procedure  of  choosing  the  most 
probable  configuration  to  represent  the  system. 


Let  us  therefore  begin  by  asking  what  is  the  most  probable  con¬ 
figuration  resulting  from  the  addition  of  many  angular  momenta,  corres¬ 
ponding  to  the  situation  in  which  many  electrons  contribute  to  the  total 
angular  momentum. 

Consider  the  N  electrons  each  carrying  an  angular  momentum  i  . 

The  resultant  angular  momentum  of  these  electrons  lies  between  zero  and  . 
If  we  consider  each  possible  arrangement  of  these  t.lectrons  equally  probable 
then  the  resultant  angular  momentum  is  shov/n  in  an  appendix  to  be  a  random 
variable  that  obeys  the  central  limit  theorem.  The  distribution  of  the 
resultant  angular  momentum  tends  therefore  toward  a  Gaussian  distribution 
centered  at  zero  angular  momentum.  This  hovjever  is  true  for  all  and 
since  the  distribution  of  a  finite  number  of  random  variables  each  of 
which  is  Gaussian  distributed  is  itself  Gaussian,  we  have  the  result  that 
the  distribution  of  resultant  angular  momentum  is  peaked  at  zero  angular 
momentum. 

In  order  to  estimate  how  the  sharpness  of  the  distribution  depends 
on  the  number  of  electrons  we  note  that  the  various  possible  vector  ad¬ 
ditions  of  individual  angular  momenta  are  equivalent  to  the  various  possible 
trajectories  in  a  random  walk  problem.  Chandrasekhar^^  gives  as  the 

probability  distribution  of  the  sum  of  a  large  number  N  individual 

•  o 

random  steps  of  mean  square  size 
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The  variance  of  the  distribution  is  . 


III.  23 


.  /,  T 
0  **  /  N  *  Ji 


while  the  maximum  possible  angular  momentum  is 

III.  24  L  =  N  /'i 

*  max  / 


V/e  see  therefore  that  the  distribution  of  the 
momentum  becomes  more  sharply  peaked  about  the  origin, 
momentum  with  increasing  number  of  particles  since 

.-1/2 


III.  25 


=:  N 


max 


resultant  angular 
i.e.,  zero  angular 


We  shall  therefore,  in  attempting  to  describe  the  situation  in 
which  many  electrons  contribute  to  the  angular  momentum  consider  only  the 
spherically  symmetric  term  in  III.  15  and  in  the  sum  on  the  right  hand  side 
of  III.  21. 

In  that  case  III.  21  becomes 


III-  26  I  “"'c 


47Td 


where 


2  2  9  r  ^ 

III.  27  KoPTTi(Rpa))=faJ-pp-il(Ul)t  +  Ze  -  2[v(R-r)n  /  (r)dr] 

_ r p  -  000  —  — '  Q 

2ni  2mR2  [r] 

'i'  .  (R) 


(R) 
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The  zeroth  order  density  n  (R.p)  is  related  to  G  (R,?) 

ooo  — ooo  — 


by 


"I-  28a  n™(Rp)  =  -i/  ^  a  8o„o<liP“> 


and 


III,  28b 


"ooo^l)  =  /  (^)3  ^oo^i’P)  Vo 


Vie  have  therefore, 


-  ,  .  n  ^47t  |d  I 

III.  29  R,paj)  =  --  --^TT 

ooo  aj-E(Rp)l) 


to  be  integrated  in  the  complex  to  plane.  For  the  case  at  hand,  i,e,,  an 
atom  in  its  ground  state  the  path  of  integration  lies  just  below  the  real  o) 
axis  from  -oo<Q)<p  ,  crosses  the  axis  at  to=vi  and  lies  just  above  the 
real  axis  from  y<oj<“  •  The  constant  p  is  the  chemical  potential.  The 
contour  is  completed  by  the  upper  semicircle  at  infinity.  The  derivation  of 
the  contour  is  discussed  fully  in  an  appendix  to  the  Baraff  and  Borowitz 
paper. 

The  result  of  the  integration  in  the  oj  plane  is 

-3 

III.  30a  n  (R.p)  =  4iTpn  E(RpJ!.)<p 

III.  30b  =  0  E(Rpil)<u 

The  spatial  density  may  be  obtained  from  n  (R*p)  by  means  of 

OOO  “ 


2h 


equation  III#  28b  which  yields 


III#  31 


n  (R)  = 
ooo  — 


1 


(  2m(iJ-<ji)'-  - 


JI()1+1)r2 

r2 


3/2 


) 


where 


-Ze^  - 

III#  32  <()(R)  =  — ^  +2  /v(R-r)n  (r)dr 

Equation  III#  31  in  conjunction  with  Poisson’s  equation  has  the 
structure  of  the  Sessler-Foley  modification  of  the  Thomas-Fermi  model. 

This  is  due  to  the  appearance  of  centrifugal  potential  on  the  right  hand 
side*  One  may  however  a.sk  why  such  a  term  appears  since  the  Green’s 
function  used  in  the  derivation  was  taken  to  be  spherically  symmetric. 

It  is  present  because  of  the  angular  dependence  of  the  transform  of  the 
kernel.  The  question  remains  however,  as  to  what  value  of  i  is  to  be  used  in 
III#  31#  To  bring  the  model  into  coincidence  with  the  Sessler-Foley  model 
we  shall  choose  i  such  -that 

J(J+l)n2  =  / 

where  J(J+l)h^  is  the  square  of  the  angular  momentum  of  the  atom.  Thus 
has  the  significance  of  an  angular  momentum  per  particle. 

As  pointed  out  by  Sessler  and  Foley,  this  model  corresponds  to 
spherical  electron  density  executing  a  rigid  rotation.  The  model  does  not 
yield  any  of  the  characteristic  properties  of  the  periodic  table  and  the 
essentially  classical  fashion  in  which  the  angular  momentum  appears  is 
exhibited  by  the  vanishing  of  the  quantum  mechanical  expectation  value  of 


the  angular  momentum. 
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In  the  Hartree-Fock  approximation  the  quantum  mechanical  ex¬ 
pectation  value  of  the  angular  momentum  is  given  by 


III. 


33  -i/LimL^CDG  (11' )d(l)t 
1'->1  1 


l2ii_/Lim2L(l)-L(2)rG,  (11'  )G,  (22'  )-G^  (21'  )G^  ( 12 '  )  ]cj(l )d(2  ) 
2  i'->-l  ^1  1  1  1 

2 '-►2 


which  may  be  evaluated  by  the  insertion  of  the  form  of  the  one  particle 
Green’s  function  employed  in  our  derivation  of  the  Sessler-Foley  result, 

i.e,  , 


Ili!il(t  -t  ) 

III.  34  G,  (r,t,r„t„)=, — (Rpw)?"  (R  )e  n  12  dpdaj 

“  —  V  no  —  poo 


1-1  1-2  2'  (2TTli)^^  00 


If  III.  34  is  inserted  in  III.  33  the  expectation  value  of  the  angular 
momentum  is  found  to  vanish  because  of  the  zero  eigenvalue  of  the  angular 


part  of  (R). 

V/e  see,  therefore,  that  if  nothing  is  known  about  the  internal 
distribution  of  the  angular  momentum  we  are  led  to  the  Sessler-Foley  result 
as  was  to  be  expected  on  the  basis  of  the  Landau  and  Lifshitz  theorem 
derived  above#  We  are,  however,  in  possession  of  a  great  deal  of  quali¬ 
tative  information  about  the  internal  distribution  of  angular  momentum  in 


an  atom.  We  would  therefore  be  led  to  expect  that  if  we  could  include  some 
of  this  information  we  would  be  led  to  a  better  model.  In  the  context 
of  the  present  fornalism  this  amounts  to  choosing  a  more  appropriate 

A 

approximation  than  Fourier-Hankel  transform  of  the  one 

particle  Green’s  function. 
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There  is  a  relatively  simple  form  of  the  one  particle  Greenes 
function  that  corresponds  physically  to  the  case  of  one,  two  or  three 
electrons  outside  a  closed  shell  and  has  the  virtue  of  yielding  the 
proper  expectation  value  of  the  angular  momentum.  A  trivial  modification 
of  the  form  of  this  Green’s  function  will  be  seen  to  correspond  physically 
to  electrons  missing  from  a  closed  shell. 

The. one  particle  Green’s  function  is  taken  to  be  of  the  form 


III.  35a  =  G°(r^t^r2t2 )  t  G^(£it^r2t2) 


where 


III.  35b  G  (r  t  r  t  )  =  i 


(r  t  r  t  )  =  ^ — :rvtr/^  (Rpo))'^  (r*  )e  n 

1  ~1  1—2  2  (2-(\'K)  ^  oo  —  poo  —1 


dpdoJ 


and 


III.  35c  G'"(r  t  r  t  )  =  ^ 


I±^(tj^-t2) 

-/g  (Rpo))'}'  (r  )e  dpdo) 

K  In  —  pLm  —1 


1  —1  1—2  2  (2iTn)^  Im'— ^  '  pLm' 

Since  the  normalization  on  the  one  uarticle  Green’s  function  is 


III.  36a  ^^^1  "  ^ 

a  suitable  normalization  on  the  quantities  just  introduced  is  given  by 


III.  36b 

o 

-i/G.j  ~ 

n=0,l,2,3 

III.  36c 

L 

•i/Gj^  (rj^t^:r^t^^)d£  =  n 
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vie  likewise  assume  that  the  two  particle  Green's  function  may 


be  decomposed  into 


Iir.  37a  G2(11':22')  ■=  G°  (ll':22')  +  G^  (ll':22') 


where 


III.  37b  G2(11':22')=  -  n/S  (Rpoo  :R'p'oou)' (.r  )'!/"  (r)e  R  ^'*^2“^2’^dpdp’dtDda) 

(2ttR)'^  2  -  -  poo  —1  p'oo  -2 


(2ttR) 


poo  —1  poo  —2 


III.  37c 


Gj  (11': 22')  = 


1  ,v  •••  ■lf^^^l"^l^“li^^^2"'‘^2^ 

( 2^mn^'LL^ ^ /S ‘■5.’  P '  Lm '  o)'  ) ) y“,Ln^  )e  dpdpUuKla/ 


The  C. .  appearing  in  III.- 37  are  Clebsch-Gordon  coefficients. 


In  choosing  appropriate  normalizations  for  the  two  terms  in 


III.  37a,  we  observe  that  the  normalization  of  the  two  particle  Green's 


function  should  be 


III.  38  jG2(ll:22)d(l)d(2)  = 


corresponding  to  the  number  of  distinct  pairs  that  one  may  choose  from 


N  objects.  V/e  shall  therefore  choose  the  normalization  of  the  two  terms 


in  the  two  particle  Green’s  function  as  follows: 


III.  39a  Jg°  (11:22 )d(l)d(2)  =  (N-n)(N-n-l) 

2  •'  2  2 


III.  39b 


/g^  (ll:22)d(l)d(2)  = 
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This  normalization  corresponds  to  the  separate  anti-symmetrization 
of  the  core  and  the  contents  of  the  outer  shell. 

If  v/e  apply  III.  33  to  the  forms  of  the  one  and  two  particle 
Green *s  functions  chosen  here  we  obtain 


III.  40  -i/LimL^(l)(G°(ll'  )+gIj'(11'  ))d(  1  )tl^~/Lim2L(l) •  L( 2 )[G° (11 ' 22  '  )+G2  (11 ' 22 '  ) )d(l )d(2 ) 
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Substituting  III.  35  and  III.  37  in  III.  40  and  remembering  that 


we  obtain 

III.  41  h2L(L+l)(-i/G{-'(ll)d(l)  J-n2L(Ltl)  2  Jg^  (1122  )d(l  )d(2 ) 

=  L(L+1)r2  ( n-n(n-l)) 

2 

V/e  thus  obtain  a  zero  angular  momentum  for  the  case  of  no 

2 

electrons  outside  a  spherical  core  and  an  angular  momentum  of  L(L+1)R 
for  the  case  of  either  one  or  two  electrons  outside  a  spherical  core. 

In  the  case  of  three  electrons  outside  a  core  the  angular  momentum  is 
seen  to  vanish  corresponding  to  the  fact  that  the  elements  of  the  principal 
groups  in  the  fifth  column  of  the  periodic  table  have  three  electrons 
outside  a  spherical  core  which  combine  to  give  an  S  state  for  the  ground 
state. 

Having  found  a  form  for  the  one  particle  Green’s  function  that 
describes  one,  two  or  three  electrons  outside  a  closed  shell  and  that 
yields  the  proper  quantum  mechanical  expectation  value  of  the  angular 
momentum  we  proceed  as  before  to  obtain  an  equation  for  the  electron 
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density*  Equation  III.  31  becomes 


III.  42 


In  this  case  the  L  appearing  in  the  Fourier-Hankel  transform  is  taken  to  be 
the  same  as  that  of  the  term  ,  As  before  we  may  evaluate  the  first 
integral  obtaining 


III.  43  Y,  (n  )K  G  =  oLm  ooo 

rb'  0  Lm  0  0  0  — - 

i  *  ‘  47TP 


^  KG  faj-E(RpL)  )  Y  ,  (R)G  (R*pu) 

nnn= '  — .  nLin  —  ooo  — 


*-(■  77  D 


8 


In  order  to  evaluate  the  second  integral,  v/e  note  the  relationship 
between  the  Fourier  transform  and  the  Hankel  transform 


in.  20b  d^(RE„)  =  ^T;;,,,(!lp)G^t„(Sp«) 


where  the  last  step  was  obtained  using  the  explicit  form  of  Using 

the  expansion  of  a  product  of  spherical  harmonics 


,  (  )m^ 0 ]_ )  % 

III.  44  Y„  (fi  )y'"'  (Q  )  =  Z  '  v;: C  (Xooo)c^  (Ay,m,-m)Y  («  )  =  Z  (Q  ) 

Jim  p  am  o  Xu'^  /il+l  ii  ll  ^4  P  Ay  Lm  Xy  p 


where  the  C  are  the  usual  Clebsch-Gordon  coefficients.  Inserting  III.  20b, 
li 

III.  44  in  III.  42  yields  for  the  second  integral  in  the  latter  equation 

Y,.  (n  )!<  =  G,  p  ?  B'^^k  ,  (Rpu) 

p  Lm  p  o  o  oLm'  Ay  Lm  oA  y  — 

=  G  .  (R)(  uj-E  ] 

U  Lm  pA  y  —  ^  ^ 


III.  45 
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where 


III.  45a 


^  ""  r  B  (r) 

Xp  Lm  pXp  - 


Thus  eq.  III.  42  becomes 


III,  46  _ 

P  ,  (R)  1^-e'  u-E(RpL) 

pLm  —  — 

Since  G  is  of  the  order  of  N  times  larger  than  G  ,  we  may 
ooo  oijTn 

solve  first  for  the  core  density  in  the  absence  of  the  outer  electrons 
obtaining  as  in  III#  29  (the  Sessler-Foley  model) 


III#  ^7 


G  (Rpoj) 
ooo  — 


H*^47TP 

co-E(R*pLT 


Since  the  functional  variation  of  particles  is  the  same 

as  for  the  N  in  the  Sessler-Foley  model  we  may  write  III#  46  as 


III.  48a 


G  ^  (Rp<i)) 
oLm 


(R) 

X  y  Lm^  pX  y  ~ 

— — m~ 

■pLm  - 


s  n 
N 


nil 


and 


III. 


48b 


G  ( Rdu ) 
ooo  — 

Integrating  as 


_  N-n  R”^4Trp 

“F'  ‘“u-E 

before  over  the  complex  energy  plane  we  obtain 


4TTpn~^n 

N 


T 

pLm  — 

'Fvy  W 

PLm  pXp 


III.  49a  n 


;  E  <  P 
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III.  49b  "ooo^-*P^  ¥  *  4ttpR-3  ;  E  <  p 


The  spatial  density  for  the  core  is  obtained  by  integrating  over 
momentum  yielding 

3/2 


Pmax 

III.  50  n  (R)  =  f  (R^r  N-n  1  J  .  L(L+l)n^i 

ooo  -  ^oC2^  N  P°°  -  N  - ^  J 

to  lowest  order  in  h. 

We  may  obtain  a  spatial  density  from  III.  49a  that  corresponds 
to  what  one  might  expect  on  physical  grounds  by  approximating  the  denominator 


by  Y  (R)  »  In  that  case 
poo  — 


Pmax 


2L+3 


III.  51  =  '^'pLm’i'pLm 

2TT2(2L+3)h2L+3 


‘  OOO 


to  lowest  order  in  h.  It  is  seen  that  in  the  case  of  vanishing  angular 
momentum  the  present  case  reduces  to  the  Thomas  Fermi  model. 

In  the  event  that  the  atom  may  be  approximately  described  by 
one,  two  or  three  electrons  outside  a  closed  core  carrying  the  entire 

angular  momentum  we  have  seen  that  it  is  possible  to  construct  a  form  for 

the  one  particle  Green’s  function  that  in  the  Hartree-Fock  approximation 
leads  to  the  proper  quantum  mechanical  expectation  value  of  the  angular 
momentum.  The  equations  for  the  spatial  electron  density  that  are 

recovered  from  this  form  for  the  Green’s  function  differ  from  those  obtained 


by  Sessler  and  Foley. 
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We  may  treat  the  cases  of  one,  two  or  three  electrons  missing  from 
a  closed  shell  by  slightly  modifying  our  previous  procedure*  We  let  n  -►  -n 
in  III*  36  and  III*  39  thereby  obtaining  for  the  expectation  value  of 
the  angular  momentum  in  III*  41 

III.  52  L(L+1)r2  (_n  _ 

Then  once  again  for  no  electron  missing  from  a  spherical  core 
we  obtain  a  vanishing  angular  momentum.  For  the  case  of  either  one  or  two 
electrons  missing  from  a  spherical  core  we  obtain  the  proper  angular  momen¬ 
tum  with  n  =  -1  and  n  =  -2  respectively*  In  the  case  of  three  electrons 
missing  from  a  closed  core  we  return  once  again  to  column  five  of  the 
periodic  table  with  a  ground  S  state* 
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IV.  QUADRUPOLE  COUPLING  CONSTANT 

The  present  model  has  been  derived  not  in  an  attempt  to  calculate 
detailed  atomic  properties  but  rather  to  illustrate  that  a  statistical 
model  including  available  information  about  the  internal  behavior  of  the 
system  will  yield  better  results.  V/e  have  therefore  chosen  to  compare 
a  calculation  based  on  the  present  model  with  a  similar  one  done  by  Sessler 
and  Foley 

The  quantity  to  be  calculated  is  the  quadrupole  coupling  constant 
IV.  1  q  =  -/  n(R)dR 

r3  - 

The  importance  of  this  quantity  lies  in  the  fact  that  the  quadrupole 
moment  of  the  nucleus  is  not  measured  directly  but  rather  through  a  quadrupole 
splitting  in  atomic  spectra.  The  relationship  between  the  quadrupole  coupling 
constant  and  the  measured  quadrupole  splitting  is 

IV.  2  hb  =  a^qQ 

where  Q  is  the  nuclear  quadrupole  moment  and  b  is  the  measured  frequency 
of  the  splitting. 

For  the  case  of  the  quadrupole  coupling  constant  of  the  halogens, 
we  may  consider  the  core  density  to  a  first  approximation  to  be  spherically 
symmetric,  the  only  contribution  to  the  quadrupole  coupling  constant  coming 
from  the  missing  electron.  Since  the  sign  of  the  second  term  in  III.  35a 
will  differ  in  the  case  of  a  single  electron  missing  from  a  closed  shell 
as  compared  to  a  single  valence  electron,  we  see  that  this  model  is  capable 
of  accounting  for  a  sign  change  in  the  quadrupole  coupling  constant  at  a 
closed  shell  whereas  the  Sessler-Foley  model  was  not. 
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Neglecting  the  centrifugal  potential  in  III*  51  we  obtain  for 
the  case  of  the  halogens 


<3  cos  0-l>r)  , 

IV.  3  q  =  •  ^3/2 


2ii2n2  L+3(2l+3)z 


2m4, 


2L+3 
2 


dR 


min 


The  expectation  value  of  3cos^0-l  in  the  ground  ^3^2 
halogens  is  -2/5*  Noting  that  for  the  present  case  t=l  and  transforming 
to  universal  Thomas-Fermi  coordinates  we  obtain 


IV.  4 


q  = 


-4/2 


25tt" 


X  . 

min 

In  order  to  arrive  at  a  reasonable  lower  limit  for  the  integration 

we  note  that  the  lowest  limit  of  the  validity  of  the  Thomas-Fermi  model 
3.6 

is  a^/Z  ,  which  is  the  radius  of  the  first  Bohr  orbit  in  the  field  of  a 
nucleus  of  charge  Ze.  That  very  little  is  lost  by  cutting  off  the  integral 
at  this  point  can  be  seen  by  considering  what  fraction  of  the  charge  of  an 
electron  in  the  lowest  p  state  in  hydrogen  lies  within  a  sphere  of  radius 
The  charge  inside  a  sphere  of  radius  in  this  case  is  approximately  ,004e 
which  is  larger  than  the  amount  of  charge  found  within  the  first  Bohr 
orbit  for  higher  p  states  and  for  higher  angular  momentum  states  as  well. 

In  addition,  for  higher  Z  this  fraction  decreases  as 

Using  Sq/Z  as  a  lower  limit  we  obtain  for  the  quadrupole  coupling 
constants  of  the  halogens 


q  a 
^Dr  o 


q  a  ^ 
I  o 


3  - 


-2.4 

-21. 

-51. 
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These  figures  are  within  an  order  of  magnitude  of  the  accepted 
values,  the  value  for  chlorine  being  in  good  agreement  with  the  accepted 
value  and  that  for  iodine  six  times  too  large.  The  value  of  |q|  however 
is  seen  to  increase  with  Z  in  contrast  to  the  Sessler-Foley  model  where  it 
varies  as 

The  primary  difficulty  in  calculating  the  quadrupole  coupling 
constant  in  any  kind  of  a  Thomas>“Fermi  model  is  that  the  density  is  related 
to  the  self-consistent  potential  and  not  to  the  square  of  a  one-electron 
wave  function.  In  the  vicinity  of  the  origin  the  square  of  the  one-electron 
wave  function  in  all  but  S  states  is  small,  in  contrast  to  the  self- 
consistent  potential  which  becomes  increasingly  Coulombic  in  character  as 
one  approaches  the  origin.  However,  it  is  precisely  this  region  near 
the  origin  that  is  weighted  most  heavily  by  the  quadrupole  coupling 


constant. 
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V.  SUMMARY  AND  CONCLUSION 


We  have  shown  how  the  inclusion  of  additional  information  in  a 
statistical  model  can  lead  to  more  desirable  results.  In  particular 
in  the  case  of  the  atom  with  angular  momentum  the  specification  of  the 
internal  structure  leads  to  a  model  that  differs  from  the  Sessler-Foley 
result# 

We  believe  that  the  present  model  offers  a  better  starting  form 
for  the  density  in  detailed  self-consistent  machine  calculations  of 
atomic  structure#  In  addition,  the  model  should  prove  useful  in  investigat¬ 
ing  at  least  the  qualitative  properties  of  the  atom  that  vary  in  the 
characteristic  way  of  the  periodic  table. 

It  is  also  anticipated  that  the  approach  used  in  the  present  paper 
may  be  utilized  fruitfully  in  any  self-consistent  statistical  model  of 
the  nucleus. 

More  generally,  however,  we  feel  that  any  statistical  treatment 
of  a  physical  system  may  be  improved  if  a  suitable  formalism  is  found 
in  which  empirical  information  about  the  system  may  be  inserted. 
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APPENDIX 

THE  CENTRAL  LIMIT  THEOREM  AS  APPLIED  TO  THE  ADDITION  OF  MANY  ANGULAR  MOMENTA 


The  various  possible  values  of  the  resultant  angular  momentum  of  N^ 
electrons  each  carrying  an  angular  momentum  Z  are  given  by  the  absolute 
value  of  the  algebraic  sum  of  N^  terms  each  of  which  may  take  on  the  23.+1 
integral  values  between  £  and  -il#  In  order  to  determine  the  distribution 
of  the  resultant  angular  momentum  we  make  use  of  the  central  limit  theorem 
which  states  that  the  sum  of  n  random  variables  of  mean  zero  tends  to  a 
Gaussian  distribution  provided 

a)  the  random  variables  possess  absolute  moments  of 
order  2  +  6  >2 


and 


b)  the  Quotient 


(i) 


2+6 


oj  =  i=l 
n 


1+6/2 


where  is  the  (2+6)^^  moment  of  the  i"^^  variable  and  B_  is  the  mean 

2+6  “ 

square  fluctuation  of  the  sum  of  the  first  n  variables,  i.e*,  the  mean  square 


fluctuation  of  the  quantity  E  k^ 

i=l 

For  a  fixed  i  the  distribution  function  is 


f(k)  = 


2i^l  I  itik 
m=-£ 


The  third  absolute  moment  exists  since 


i 

1  E 

2  ^+1  m=-it 


k 

3|  .  2 


k,m2^+l 


^mk 


m 


2t+l  Tn=0 


m 


_ 2_  (Uifi)) 

2U1  2  ^ 


The  mean  square  fluctuation  of  is 

o  o’  2 

kT  =  ^  k.f(k,)  k.  £^=I  2^=0 

1  k  1  1  k,m  i2i^tl  .  2ifl 


m=-il 

^  2  ,  2^/,  V  ,2  6Tnki  t 

ki  =  ^  =  Z  k^ 

k  k,Tn  m=-Jc 

£ 


Thus  Ak.^  =  liilil  and 
i  3 


j  .2  =  »(yf) 

m=0  ^ 


=  nAk^ 


2  -  nJl(£.+l) 


The  quotient  therefore  becomes 

r?.(£+i) 


w_  = 


^  22,+!  n/o  0/9 


which  goes  to  zero  as  with  increasing  n. 

Since  the  distribution  of  angular  momentum  satisfies  the  central 
limit  theorem  it  tends  in  the  limit  of  large  n  to  a  Gaussian  distribution 


which  was  to  be  shown. 
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